Radiation reaction in quantum field theory 
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We investigate radiation-reaction effects for a charged scaiar particie acceierated by an externai 
potentiai reaiized as a space-dependent mass term in quantum eiectrodynamics. fn particuiar, we 
calculate the position shift of the final-state wave packet of the charged particle due to radiation at 
lowest order in the fine structure constant a and in the small h approximation. We show that it 
disagrees with the result obtained using the Lorentz-Dirac formula for the radiation-reaction force, 
and that it agrees with the classical theory if one assumes that the particle loses its energy to radi- 
ation at each moment of time according to the Larmor formula in the static frame of the potential. 
However, the discrepancy is much smaller than the Compton wavelength of the particle. We also 
point out that the electromagnetic correction to the potential has no classical limit. (Correction. 
Surface terms were erroneously discarded to arrive at Eq. 1591 1. By correcting this error we find 
that the position shift according to the Lorentz-Dirac theory obtained from Eq. Ill 2t is reproduced 
by quantum field theory in the h — > limit. We also find that the small V(z) approximation is 
unnecessary for this agreement. See Sec. IVllI ) 
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I. INTRODUCTION 



Although classical electrodynamics is a well established 
theory, there has been some controversy over the nature 
of the reaction of a point charge to its own radiation. 
The standard equation for the radiation-reaction force, 
the Lorentz-Dirac equation P, 0] , admits runaway solu- 
tions which describe the charged particle accelerating on 
its own and reaching the speed of light in a very short 
time. (See, e.g., Refs. 0,0 for a review of the Lorentz- 
Dirac theory.) Although it has been shown that there 
is a unique solution that does not exhibit the runaway 
behaviour for any bounded time-dependent force acting 
only for a finite interval of time [f| (see also Refs. 0,0), 
this solution violates causality. It is generally believed 
that these problems will be absent once the finite-size 
effect is properly taken into account. For example, an 
analysis of a charged sphere has led Yaghjian 7] to the 
conclusion that the radiation-reaction 4-force is modified 
at the onset of acceleration in such a way that there is no 
causality violation in the unique non-runaway solution. 
There is also an approach based on "reduction of order" 
[H which excludes runaway solutions without introducing 
causality violation. In this method the radiation-reaction 
force is re-expressed in terms of the external force. This 
approach has recently been advocated by several authors 

mm 

Thus, it may be fair to say that the aspects of the 
Lorentz-Dirac theory that were traditionally regarded 
as problematic have been clarified. However, there are 
other counter-intuitive aspects in the Lorentz-Dirac the- 
ory which are not discussed very often. These aspects 
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persist in reformulations of the Lorentz-Dirac equation 
mentioned before. For example, the radiation-reaction 
force vanishes if the acceleration of the charge is constant. 
Although this is not so serious a problem as existence of 
runaway solutions or causality violation, it is certainly 
counter-intuitive since the charge radiates energy away 
continuously while being accelerated. A related fact is 
that the radiation-reaction force makes the kinetic energy 
of the charge increase if the acceleration of the charge in- 
creases in time, as we will see in the next section, despite 
the fact that the charge is radiating in the meantime. 

In view of these counter-intuitive features of the 
Lorentz-Dirac theory, it is natural to ask whether or 
not this theory reproduces the classical limit for a point 
charge in quantum electrodynamics (QED). (Moniz and 
Sharp have studied QED in this context and concluded 
that the Lorentz-Dirac theory is reproduced in the clas- 
sical limit by studying the Heisenberg equations for an 
extended charge and taking the zero-size limit 11].) It is 
well known that the loopwise expansion in quantum field 
theory is an expansion in powers of H |12l |l3( . Therefore 
the classical limit h — > in this sense is the tree approx- 
imation. However, one should not use this limit in the 
problem at hand. For example, this limit is taken with 
the mass parameter mc/h held fixed (see, e.g., Ref. 14|). 
Thus, the mass m also tends to zero in this limit. In 
the radiation-reaction problem the mass as well as the 
electric charge e must be kept fixed in the limit h — > 0. 
Therefore, the limit we would need to consider is cer- 
tainly not the tree approximation and, in fact, does not 
exist as emphasized recently by Julia |15| . This fact can 
readily be seen by noting that the fine structure constant 
a = e? /Antic 1/137 is of order ftr 1 . Thus, perturba- 
tion theory would break down in the h — ► limit with e 
fixed. For example, the one-loop correction to the elec- 
tron magnetic moment, a/2-K, would diverge despite the 
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fact that it is a very small fixed number. 

Although the limit K — ► with m and e fixed would 
be divergent, one can still compare some quantities in 
QED with classical ones at order a. At this order in a, 
physical quantities of interest are of finite order in hr 1 , 
and there are some quantities of order H , which can be 
compared with the corresponding quantities in classical 
electrodynamics. We note, however, that these quantities 
are much smaller in reality than some quantities which 
vanish in the limit H — > 0. The following example illus- 
trates this paradoxical feature. Let w be the frequency 
of an electromagnetic wave. Then the energy of a single 
photon htu tends to zero as H — > 0, but ahuj — e 2 uj/4irc 
remains finite. Thus, the quantity afkv is meaningful in 
classical electrodynamics even though it is much smaller 
than Hlo. 

In this paper we examine the behaviour of a wave 
packet of a charged scalar field undergoing acceleration 
due to an external potential realized as a position de- 
pendent mass term to first order in a. In particular, 
we calculate the change in the position of the final-state 
wave packet, the position shift, as a result of radiation 
using the WKB approximation. We demonstrate that 
the h — > limit of this quantity disagrees with the cor- 
responding result in the Lorentz-Dirac theory. Then we 
show that it agrees instead with the result in classical 
electrodynamics obtained by assuming that the kinetic 
energy of the particle is lost to the radiation given by 
the classical Larmor formula at each moment of time in 
the rest frame of the potential. We also point out that 
the electromagnetic correction to the potential is of or- 
der ft -1 and, therefore, cannot be calculated in classical 
electrodynamics. The rest of the paper is organized as 
follows. In Sec. |H] we review the Lorentz-Dirac theory 
and calculate the position shift of a charged particle due 
to radiation in the nonrelativistic approximation in this 
theory. In Sec. IIHI we calculate the transition amplitude 
for a charged scalar particle emitting a photon using the 
WKB approximation. This is used in Sec. IIVI to calcu- 
late the position shift in scalar quantum electrodynam- 
ics. Here, we closely follow an unpublished paper |l6j| . 
Then in Sec. E] we examine the one- loop effect (with no 
emission of photons) which results in a correction to the 
external potential. We summarize our results and make 
some concluding remarks in Sec. IVII The calculations 
are performed using conveniently chosen physical polar- 
ization vectors for photons, but they are shown to agree 
with the results obtained without choosing physical po- 
larizations in Appendix A. A formula used to compute 
the position shift is justified in Appendix B. Our metric 
signature is H . 



II. THE POSITION SHIFT IN THE 
LORENTZ-DIRAC THEORY 



rally arises in classical electrodynamics. (See, e.g. Ref. 
.4] for a derivation.) Suppose that a classical charged 
particle with charge e and 4-velocity u^ L is accelerated 
by an external 4-force F^ xt . [The 4-velocity of a parti- 
cle moving in the z-direction with speed v is (7c, 0, 0, 71*) 



with 7 = (1 - v 2 /c 2 )- 1 / 2 . Note that u^u h 
4-velocity satisfies the following equation: 
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where m is the rest mass of the particle and where r is the 
proper time along the world line of the particle. The par- 
ticle emits electromagnetic radiation with 4-momentum 
P 1 * = (E/c, P), where E is the energy and P is the mo- 
mentum, given by the relativistic generalization of the 
Larmor formula: 



dP^ 1 

~dV 
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(2) 



where a = ^—ii^u^ (with ii^ = du^/dr) is the proper 
acceleration of the particle (see, e.g., Ref. 0,3). Since 
it is the particle that radiates, one expects that its 4- 
momentum should be reduced by the amount carried 
away by the radiation. It is natural to describe this ef- 
fect in terms of a 4-force. Thus, one is led to modify the 
equation of motion of the particle as 



du^ „ 



(3) 



where is the radiation-reaction 4-force which repre- 
sents the back reaction of the radiation on the charged 
particle. It may seem that the 4-force should be equal 
to the negative of the right-hand side of Eq. J5J because 
of energy-momentum conservation. However, this can- 
not be the case because this 4-force would not satisfy the 
condition u^K^ = 0, which is a consequence of the equa- 
tion u^ii^ = 0. Therefore, one is led to require only that 
the total loss of 4-momentum of the particle be equal 
to the total 4-momentum radiated away assuming that 
the acceleration occurs only for a finite time. Thus, one 
requires that 



6ttc 



a 2 u^ + 



d(y_ 

~dT 



(4) 



for some vector because the second term does not 
contribute to the total change in the 4-momcntum, 
f^K^dr. By using the condition u^K 11 = one ar- 
rives at the following equation: 



dC» 



dr 67rc 3 
Noting that a 2 — —u^u^ — u^il^, one finds 



(5) 



We first motivate the standard Lorentz-Dirac theory 
briefly. We do not derive it but describe how it natu- 



dC 
~dV 
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67rc 3 dr 2 



dr 
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where dC^/dr is orthogonal to u M . Letting C 11 = 0, one 
arrives at 



67rc J 



d 2 u^ 
~dP r 



1 



(7) 



This is known as the Abraham 4-vector [13, and Eq. © 
with this 4-force is the Lorentz-Dirac equation. 

To illustrate counter-intuitive features of this equation 
let us consider one-dimensional motion parametrized by 
the rapidity f3 with u M = (it , it 1 ) = (ccosh/3, csinh/3) 
and F e ^ t = (F ext sinh/3, # oxt cosh/3). Then Eq. © reads 



-J— — -Text + 

ar mc 



d 2 /3 



6irmc 3 dr 2 



(8) 



The second term represents the radiation-reaction force. 
Notice that it vanishes if the acceleration is constant 
and is positive if the acceleration increases. Thus, the 
radiation-reaction force pushes the particle forward if the 
acceleration increases even though the particle is radiat- 
ing. This is rather counter-intuitive. 

Now, let us describe the classical counterpart of our 
model. In this model a particle with mass m and charge 
e moves in the z-direction from z ~ — oo to z = +oo and 
is accelerated by a static potential V(z) whose deriva- 
tive is nonzero only in a finite interval. We define 
V-oo = V(— oo) and let V(+oo) = 0. We work in the 
nonrclativistic approximation. Thus, the energy is given 
by \mv 2 + V(z), where v is the velocity of the particle. 
We assume that the kinetic energy is much larger than 
the potential energy and we retain only the terms up 
to second order in V(z). Note that the particle moves 
in the positive z-direction forever under this assumption. 
We assume also that at t = the particle has passed 
the region with V'(z) ^ and that it moves at con- 
stant velocity for t > 0. Let us denote the position and 
velocity of the particle without radiation by z and v, re- 
spectively, and define Sz = z — z and Sv = v — v. One 
has 5z\ t= -oo = and £u|t=-oo = by definition. We 
will calculate the position shift due to radiation at t = 
denoted by 5z\t=o using the Lorentz-Dirac equation. We 
emphasize that this position shift is measured after the 
acceleration because there is no acceleration for t > by 
assumption. 

The rate of change in the energy of the particle accord- 
ing to the Lorentz-Dirac theory is 



d_ 

dt 



—mv 2 + V(z) 
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vv , 



(9) 



where the dot indicates differentiation with respect to 
time t. By writing vv as d(vv)/dt — v 2 and integrating 
with respect to t, we find 



lmv 2 + V(z) = E + ^ 



vv 



6ttc 3 



v 2 dt, (10) 



where E is the initial energy. By subtracting the equation 
\mv 2 + V(z) — E from Eq. (|1(J|I we find to lowest order 



m e 



mvSv + V'(z)Sz 



e 2 „dv e 2 rt 
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d 2 dt, (11) 



where d = dv/dt. By substituting V'(z) — —mdv/dt in 
Eq. |(TJ| and dividing by v 2 , we obtain 



,2 rt 



dt V v J 67Tuc 3 dt 6nv 2 c 3 



d 2 dt. (12) 



By integrating this formula from t = — oo to £ = we 
have, to second order in V{z), 



1 9 
Sz\ t =o = gU/ log 3 I , 



(13) 



where V{ and Vt are the initial and final velocities, re- 
spectively, without radiation, and where 

1= [ If [d(t')] 2 dt'\dt. (14) 

We have replaced mv by the final momentum p = mv t in 
the second term in Eq. l|13f) because we retain the terms 
only up to second order in V(z) (and consequently to 
second order in the acceleration 5). Integrating by parts 
and using the assumption that a(t) — for t > 0, we find 



/+oo 
t[d(t)] 2 dt. (15) 
-oo 

Now, the power of radiation can be found from Eq. © 

as 



(16) 



in the nonrelativistic approximation. (This is the Larmor 
formula.) Then the position shift can be expressed as 



Sz\ t= o 



6nmc 3 



v f i r + °° 

v f log^- + - / tPJt)dt. (17) 

Vi P J-oo 



In terms of the Fourier transform a p (ui) of a(t) defined 



by 



we have 



-Hex; 



Op(u) = / dtd(t)e luJt , 



T ■ [ + °° ^ - ( \* d - f ^ 

1 = l J_ 00 ^ apH ' 



(18) 



(19) 



where A* denotes the complex conjugate of A. (The sub- 
script "p" in a p (u>) has been inserted to emphasize that 
this quantity depends on the momentum of the particle. 
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This notation will turn out to be useful later.) Thus, to 
second order in V{z) we have 



Sz\ t =o 



67ttoc 3 ^ Vj 



67175c 3 



/_r£ w ^ a ^- (2o) 



The first term is first-order in the potential because v/ — 
Vi is. ft is rather puzzling because the position shift is an 
effect caused by radiation which is second-order in V{z). 
On the other hand, the second term is clearly second- 
order. Note also that the first term would be absent if 
Vf = Vi, i.e. if V-oo = 0. In Sec. II VI we will derive the 
position shift Sz\t=o in the WKB approximation in scalar 
QED. We will find that the first term in Eq. l(2U|l is absent 
in the classical limit even if vt ^ Uj. 



where = (x,y), p± = (p x ,p v ) and 



Po/c— J m 2 c 2 + p 2 + p\ 



(26) 



with p > 0. Thus, the quantity p is the z-component of 
the momentum in the region with large and positive z for 
which the potential V(z) vanishes. The function 4> p (z) 
satisfies 



2m dz 2 



V(z) 



(j)p(z) 



P 

2m 



(j)p(z) 



(27) 



This function can be given in the WKB approximation 
(see, e.g., Ref. [5) as 



(p p (z) 



k p {z) 



exp 



K p (z')dz' + g(z) 



(28) 



III. THE TRANSITION AMPLITUDE 

We consider a complex scalar field ip(t, x) which is cou- 
pled to the electromagnetic field A^(t, x) and accelerated 
by an external potential V(z) with the properties de- 
scribed in the previous section. This model is given by 
the following Lagrangian density: 

C = [<9 M - i(e/hc)A^ ■ [d^ + i{e/hc)A^]i! 



with 



n 2 



TO c 2 + 2mV(z)}ip^ 



-F^, (21) 



with do = c~ 1 dt and = d^A^ — d^A^. The mode 
functions <f>(i,x) for the free scalar field (with e = 0) 
satisfy 

[d^ + m 2 c 2 /h 2 + 2mV{z) /fr 2 ] x) = . (22) 

The energy Pq 1 and the momentum p cl of the correspond- 
ing classical particle satisfy 



iPo/c) + (P ) + "*V + 2mV(z) = 



(23) 



(Note that the time component of the 4-momentum p c ^ is 
not pg 1 but PoV c -) I n the nonrelativistic approximation 
with m 2 c 2 > (p cl ) 2 , 2m\V(z)\, we have 



p c l m mc z + (p c1 ) 2 /2to + V(z) 



(24) 



Thus, the particle moves under the influence of the 
nonrelativistic potential V(z), and the quantum model 
given by Eq. (|21|) corresponds to the classical one an- 
alyzed in the previous section. Later we will assume 
to 2 c 2 ^ (p cl ) 2 ^ 2m|y(z)| - this condition implies 
that the particle is nonrelativistic but the kinetic energy 
is much larger than the potential energy — but we do 
not use this assumption for the moment. 

The solutions of Eq. (1221 which are relevant here can 
be written as 



(i,x) = p (z)exp 



■z(~Pot + ip± ■ xj_) 



(25) 



k p {z) = [p 2 -2mV{z)] 1 ' 2 , 



(29) 



where g(z) is the correction term of order h. We let 
g(z) = to work to lowest nontrivial order in h. The 
wave function (12811 is normalized so that 



dz 4>pi (z)* 4> v {z) = 2i:h5(p — p') 



(30) 



Strictly speaking, this formula is not correct because the 
WKB approximation is not valid for some values of p. 
However, since we use only the modes for which the WKB 
approximation is valid, the final result will not be affected 
even if we formally use this approximation for all range 
of p, as we do here. Then, the field ip can be expanded 
using these modes as 

=t*J 2^(2^)3 [^(P) $ P^) + Bt (P) $ pWl . 

(31) 

where x = (t,x), p = (p, pj_). The modes $ p (i,x) are 
orthonormal: 



ih J d 3 x$ p (i,x)* d t $ p /(i,x) = 2p (2 7 r? l ) 3 ( 5 3 (p-p'), 

(32) 

where dt = dt — dt- By imposing the usual canonical 
commutation relations on ^(i,x), 

|V>(t,x),fl^ t (t,x')] = i/ic 2 <5 3 (x-x'), (33a) 
[ip1(t,x.),dti>(t,x')] = ihc 2 5 3 (x-x'), (33b) 

with all other equal-time commutators vanishing, one 
finds 

[A(p),A\p')} = 2 Pa (2Trh) 3 6 3 (p-p'), (34a) 
[B(p),ijt(p') = 2p (2 7 rft) 3 ( 5 3 (p - p') . (34b) 
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-iwt+ik-x 



+ h.c] , (35) 



All other commutators of the creation and annihilation 
operators vanish. The electromagnetic field in the Feyn- 
man gauge can be expanded in the usual manner as 

where u> = c||k||. The operators &u(k) satisfy 

[6 M (k),6,(k')] - 0, (36a) 
[&„(k),bt(k')] = - 9Al ,2^(27r) 3 ( 5 3 (k-k'), (36b) 

where g^ v is the metric of the Minkowski spacetime. No- 
tice that we have used the wave number k instead of the 
momentum to label the modes. This is more convenient 
because the frequency and wave number of the emitted 
photon are classically well defined, the former being re- 
lated to how rapidly the acceleration changes. Thus, the 
momentum of the photon is of order h. 

We perform our calculations with physical polariza- 
tion vectors eW^k), j = 1,2, of the photon with wave 
number k, which are real and with vanishing time com- 
ponents. We choose them so that e^ z (k) — 0. It is 
shown in Appendix A that the calculations using the 
Fock space with indefinite metric in the Feynman gauge 
yield the same results. We define the transition ampli- 
tude from the initial state consisting of a charged particle 
with momentum p = (p, p± ) to the final state consisting 
of a charged particle with momentum P = (P, Pj_) and 
a photon with wave number k and polarization vector 



(k) by 



-A(j,k,P,p) 



I f d 4 x(0\b^Ck)A(P)C I (x)A^(p)\0) 



where 

6«(k) = e^>(k)^(k), 
with the following interaction Lagrangian density: 



IE 

he 



(37) 
(38) 

(39) 



leading order in h. Since we know on physical grounds 
that the energy and momentum of the photon emitted 
are of order h and, therefore, much smaller than those of 
the charged particle, we can assume that P ^> ||Pj_|| for 
the final state because of the assumption that p ^> ||p_i_|| 
for the initial state. Thus, we neglect the term propor- 
tional to (p ± + P_l) • (k) in Eq. gQJ. Then, to lowest 
nontrivial order in h we find 



We only need the transition amplitude with the initial 
state satisfying p ^> ||p_i_|| because we will use a wave- 
packet state with this condition satisfied. By a straight- 
forward calculation we find 

_4(j,k,P,p) = ec 2 ydz{(p ± +P ± ). e ^(k)0 P (z)*^(z) 

x(2Trh) 3 s 2 (P ± + hk ± -p ± ) 

xS(P + hLU - Pa ) (40) 

with k ± = (k x ,k y ), ef = (e<j) x ,eV>v) and d z =d z 
— dz- We will calculate this transition amplitude to 



-4(j, k,P,p) 



(k)G(k z ,P,p) 



x(2ttK) 3 S 2 (P 1 _ + hk ± -p x ) 
xS(P + hw-p ) 7 (41) 



where 



-ik z z 



G{k Zl P,p) = J dz<ft P (z)* [Kp(z) + Kp(z)]<f> p (z)e 

(42) 

Next we evaluate the integral G(k z , P,p) to lowest non- 
trivial order in K. We keep only the terms which are first- 
order in V(z) because G(k z ,P,p) will be squared in the 
position shift for which we need only the second-order 
terms. First we note that to first order in V(z) we can 
write 



G{k z ,P,p) 
where 



2-wPp / dzexp 



i l K(z')dz' 
o 



K(z) = [k p (z) - K P (z)} /h - k z 



(43) 



(44) 
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This follows from the fact that [k p (z)/kp(z 
[np(z) I Kp(z)] 1 / 2 — 2 is second-order in V(z). Now, for 
any function g(z), one has 



dze* 9 ^ = 



WW? 



(45) 



We use this formula with g(z) — K{z) dz + ie\z\, where 
the term ie\z\ with e > has been inserted to regularize 
the integral. The boundary terms vanish because one 
cannot have K(z) = as well as 



Pa 

P-L 



Po 
P 



(46a) 
(46b) 



which are enforced by the 8- functions in Eq. (|41|l . with- 
out violating the mass-shell conditions. Hence we have 



G{k Zl P,p) = -2i^f¥~p J 



dz 



x exp 



K'{z) 
[K(z)\ 2 

K(z')dz' 



Let us find a small h approximation for K(z). 
squaring the energy-conservation equation l|46al) , we 
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■Px 



p2 



(47) 
By 

find 

(48) 
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in the nonrelativistic approximation and to first order in 
h. By using the inequality 

|pi-Pil<||px-Px||||p± + Px|| (49) 
and Eq. (|46b|) . we find 

|P± -P±| < Hu/c- ||Px + Px|| ^.mhuj. (50) 
Using this in Eq. Ij48(l we have 

p 2 ~ P 2 ^ 2mhuj . (51) 
Hence we obtain the following approximate formula: 



K p (z) - Kp(z) 



dP 2 KP ' ~ v p (z) 



(52) 



where 



v p (z) = K p (z)/m = \J p 2 — 2mV(z)/m (53) 



is the speed of the corresponding classical particle. 
Hence, 



K(z) 



v p (z) 



k z . 



(54) 



Furthermore, since the charged particle is nonrelativistic, 
we have uj/v p {z) = c\\k\\/v p (z) ^> k z . Thus, 



K{z) 



v p (z) 



(55) 



By substituting Eq. in Eq. and letting ^JTp ; 
p, we obtain 



G(k z ,P,p) 
2ip 



dz v' (z) exp < iw I \v p {z')] 1 dz 



Let us define the function T(z) by 



Tee / [v p (z')]- l dz' +T , 
Jo 

where To is a constant. Then by using 

dz v' p {z) = dz ^ = dT 



(56) 



(57) 



(58) 



and noting that T — > ±oo as z — > ±oo because v p {z) is 
always positive, we obtain 

G(k z ,P, P ) = [ +C °dT^e^ T -^ . (59) 
w J-oo dT 

Now, Eq. (|57() implies that dz/dT = v p . Thus, the quan- 
tity T can be identified with the time coordinate of the 
classical particle if the constant T is chosen appropri- 
ately. Recall that we need to evaluate G(k z , P,p) only to 



first order in V(z). Since the acceleration dv p /dT is first- 
order in V(z), we need to find T only to zeroth order in 
V(z), i.e. we may neglect the potential. Since T = Tq for 
z = 0, the constant To is the time when the particle is at 
z = 0. Letting the particle be at z = zq for T = 0, we 
have zq/Tq — —v p = —p/m. Hence To = —mzo/p. By 
substituting this in Eq. (|59|l we find 

lip 

G(k z ,P,p) = a p (ui) exp(imojz /P) , (60) 

co 

where the Fourier transform of the acceleration, a p (u>), 
is defined by Eq. (|18|) [with the identification dv p /dT = 
a(t)]. One may replace p by P and vice versa in Eq. I|60|) 
because the difference p — P is of order h. 

It is instructive to compare this result with the tran- 
sition amplitude for radiation due to a classical current 
with the following interaction Lagrangian density: 

C Ici = (e/c)A M (i,x)^(i,x). (61) 

Here, the current density j^ix) is given by 

j°(t,x) = c5 2 (x ± )S[z-z(t)], (62a) 
j 2 (t,x) = v(t)8 2 (x±)5[z-z(t)] (62b) 

with j x (t,x) — j v (t,x) = 0, where z(t) is the position of 
the charged particle at time t and where v{t) = z'{t) is 
its velocity in the z-direction. The transition amplitude 
from the zero-photon state |0) 7 to a one-photon state 
&Wt(k)|0) 7 is 

Ai(i,k) ee i J d 4 x 7 (0|&«(k)£ /cl (x)|0} 7 

/ + oo 
dtv(t)exp[iu;t~ ik z z{t)\ . 
-OO 

(63) 

By integrating by parts with a damping factor e _e '* and 
dropping the surface terms because u> — k z v(z) ^ for 
all z, we obtain 

A cl (j, k) = ieeWfk)*, f °° d/ (<) 6XP M ~ lkzZ(t)] 



k z v(t)Y 



Since lo = c||k|| k z v(t) and tot 3> k z z(t), we find 
Ai(j,k) « W / dtv'(t)e wt 

U J-oo 

iee u *> z (k) 

UJ 



(64) 



a p (u) . 



(65) 



Notice the similarity of this expression with Eq. 160fl . 
Our result obtained in the small H approximation is 
closely related to the emission process by the correspond- 
ing classical particle. 
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by 



IV. THE POSITION SHIFT IN QUANTUM 
ELECTRODYNAMICS 

Now, let us consider an initial wave-packet state given 



r/(p)^(p)|0) 



/2p~o~(2irh) 3 ' 

The normalization condition (i \i) = 1 leads to 



d 3 p 

(27rft) : 



I/(P)| S 



1 



(66) 



(67) 



The function /(p) will be identified later as the non- 
relativistic one-particle wave function in the momentum 
representation. We choose the following form: 



/(P) = |/(p)]e" 



-ipzo/H 



with 



|/(p)| = 7V ex p 



(P-P) 2 + \\P±\ 
2(Ap) 2 



(68) 



(69) 



where N is the normalization constant. Here, Ap and 
p(> 0) are the width of the wave packet in the mo- 
mentum space and the average momentum [in the region 
with V(z) — 0], respectively. The constant zq will turn 
out to be the average z-coordinate at t = 0. We now 
assume that m 2 c 2 » p 2 >• m|V r (z)|. We also assume 
that Ap <C p. The transverse momentum is much 
smaller than p under this assumption. This wave packet 
describes a nonrelativistic charged particle moving in the 
positive z-direction with its kinetic energy much larger 
than |V(;z)|. The initial state |z) given by Eq. (|FJF)|) leads 
to the following one-photon final state according to the 
standard time-dependent perturbation theory in the in- 
teraction picture: 



\L 17) 



d 3 k 



d 3 p 



2^(2tt) 3 J 2P (2irh) 3 



d 3 p 



/2p^(27rft) 3 
x6( i)f (k)At(P)|0) 



/( P M(j,k,P,p) 



(70) 



where Pq/c = ^/m 2 c 2 + P 2 + P 2 j_. The amplitude 
A(j, k, P, p) is defined by Eq. (gTJ with G(k z ,P,p) given 
by Eq. f>UJ). By substituting Eq. |gT} in Eq. (J7UJ) and 
using, for any function F (p), 



d 3 p 

V^Po 
2pc 2 



F(p)6 2 ( Pl _ + hk ± - P ± )S(P Q + hio-p ) 



(71) 



where p on the right-hand side is expressed in terms of 
P and k through the relations l|46|l . and recalling that 



e ( 2 ) z (k) = 0, we obtain 



I/, 17) = f 



d 3 k 



d 3 P 



Hp) 



Xf 



2w(2tt) 3 J 2P (2i:h) 3 2p 
!l ^(k)G(fc,,P, P )&«t(k)At(P)|0) . (72) 



Now, from Eq. (jf> 1|> we have p — P w mhuj / P. By using 
this approximation in Eq. I|68|l we find 



f(p) ~ l/(p)l exp (~iPz /h - imujz /P) 



(73) 



By substituting this equation and Eq. I|60|l in Eq. I|72l) 
we obtain 



\f,H 



d 3 k 



d 3 P 



2w 2 (2tt) 3 J V2T^(27rfi) 3 
x |/(p)| exp (-iPz /H) e« z (k)a» 
x6( 1 't(k)A t (P)|0>, (74) 



where 



P=(p,P±) = (P + mhw/P,P, 



hk. 



(75) 



We have approximated po by Pq on the right-hand side 
of Eq. JZH. 

The emission probability can be obtained to leading 
order in h from Eq. (|74|l as 



V r = (/,1 7 |/,1 7 ) 
d 3 k 



d 3 P 

h J 2w 3 (2tt) 3 J (2irh) 3 

x|/(p)| 2 [ £ (^(k)] 2 |a p H| 2 
- 2 ' d 3 k r ... 12 



2w 3 (2tt) 



«*(k) |a F M| 2 , (76) 



where we have used the fact that the function |/(p)| is 
sharply peaked about (p,p±) = (p,0). Recalling that 
d 3 k = dtlkU} 2 dcu/c 3 , where dfik is the solid- angle element 
in the k space, and that the average of |e^ z (k)| 2 for each 
value of a; is 2/3, we find 



STthc 3 



JO 27TW 



(77) 



Notice that this probability is infrared divergent unless 
Op(0) = V-oo/p = 0. However, the expected energy of 
the photon is finite and is given by 



E r — 



37rfi,c 3 
e 2 

67TC 3 





+00 



2t:lu 



hui\dp(uj)\' 



dt [a ¥ (t)Y 



(78) 



where ap(t) is the acceleration of the corresponding clas- 
sical particle at time t. This is nothing but the Larmor 
formula (|16|l . 
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For the classical current discussed at the end of the 
previous section, the one-photon final state is 

d 3 k 



l/,l7>c 



h J 2w(2tt) 3 
e f d 3 k 



h / 2w 2 (2tt) 3 



Ai(l,k)&< 1 )t(k)|0) 7 

«*(k)a pM &«t(k)|0) 7 



The emission probability calculated from this state agrees 
with Eq. If77jl . 

Now, we derive the quantum expression of the posi- 
tion shift to order e 2 in the small h approximation and 
compare it with the result in the Lorentz-Dirac theory 
derived in Sec. |H] We assume that the constant zq in 
Eq. I|t)8|) is larger than any z with V(z) ^ and that the 
wave packet is far into the region with V(z) — at i = 0. 
First we note that the expected value of the z-coordinate 
at t = with the electromagnetic field turned off is given 
by 



\off 
/t=0 



if). 

T 



d 3 p 

(2irh) 3 



Hp) 



.9 Hp) dfjpr 

dp 



/(P) 



dp 

(80) 

where /(p) is defined by Eq. (|68|l . This formula can 
be understood as the statement that the function /(p) 
is the momentum representation of the one-particle wave 
function. The derivation of Eq. 18Ufl is given in Appendix 
B. Since 



ih 

T 

we find 



/(p)«-^/(p) 



dp 



dp 



Z . 



^o|/(p)| 2 , (81) 



(82) 



Next we calculate the expected value of the z- 
coordinate at t = with the electromagnetic field turned 
on. At lowest nontrivial order in e, the final state can be 
given as 



l/> = l/,07> + l/,l7>. 



where |/, I7) is the one-photon final state given by Eq. 
(|74|l and where the zero-photon final state |/, O7) can be 
written as 



d 3 i 



: [l + ^(p)]/(p)^(p)|0). (84) 



I/.07) 



We will not explicitly evaluate the one-loop forward- 
scattering amplitude .F(p), which is of order e 2 . The 
expected value of the z-coordinate at t = is the sum 
of the zero-photon and one-photon contributions. The 
zero-photon contribution, (z)^ , can be obtained by re- 
placing /(p) by [1 + iF(p)]/( p ) in Eq. (JHDJl. The result 



is 



v(0) 
>t=0 



d 3 p 

(2irhy 



I/p 



z (l - 2ImJT) 



dp 



(85) 



where we have written /(p) = / p and used Eq. I|81|l . The 

contribution from the one-photon state, (z)[l} , is similar 
except that one needs to trace over the photon states. 
Recall that we have calculated G{k z ,P,p) in Eq. I|60() to 
lowest order in h. We define ai(k, p) as a quantity such 
that the expression for G{k Zl P,p) in that equation will 
be valid to next order in h if e*- 1 '' z (k)aj,(w) is replaced by 



(79) it. (Introduction of ai(k, p) is necessary because (z 



and (z)[]} are of order h 1 .) Then, we have, to order h°, 



(o) 
t=o 



d 3 k 



ie 

T j 2^ 3 (2tt) 3 



d 3 P 



x 



'«' i(k P)| 2 (/p^-^/p 



(k)| 2 |/ P | 



„9a p (a;) da p {uS)* 



dP 



dP 



a p (uj) 



(86) 



where p is related to P by Eq. (|75|l . The integration vari- 
ables P can be changed to p if we change the derivative 
d/dP to d/dp at the same time. Then by using Eq. I|81|) 
again we find 



z) (1) - 
z /t=o — 



d 3 P \f P \ 2 f d 3 k 
(2nh) 3 J 2w 3 (2tt) 3 

^|a l( k,p)| 2 + ^| £ ^(k)| 2 



CL p {ui) 



„9a p (w) dd p (ujy 



dp 



dp 



dp(uj) 



(87) 



Unitarity implies that (/, 7 | /, 7 ) + (/, 1 7 | /, 1 7 ) = 1 
at order e 2 . This results in the following relation: 



2ImJ- = _ 

h 



d 3 k 



2w 3 (2tt) : 



|ai(k,p)| : 



(83) Both sides are equal to the emission probability, which is 



,(0) 



given by Eq. (|77J) at leading order in h. By adding 

and (z)!lo given by Eqs. (|H3f) and (JSSJ), respectively, using 
Eq. H88|) and recalling that the function / p is sharply 
peaked about (p, p) = (p, 0) we have 



( z )t=o 



dp 



ie 
~~2 



P=P>P_L=0 

rf 3 k | £ (l) 2(k) |2 

2w 3 (27r) 3 



d d 
d p (ui)*— a p (uj) - —d v (uj)* ■ d p (uj) 



P =p 
(89) 



Now we note that the classical acceleration depends on 
p and t through z — pt/m + Zq. [One can use the free- 
particle approximation here because the last term in Eq. 
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(|89|l . where a p (uj) appears, is already second-order in 
V(z).] Thus, 



V. SOME COMMENTS ON THE 
FORWARD-SCATTERING AMPLITUDE 



— - — 

dt dp 



Then we find 

d_ 

dp 



a p (u>) 



+oc da; 

—a p {u)e 



--a p {io) —a p (u>) . 

p p OLU 



(90) 



(91) 



From this formula and the fact that a p {— uj) 
[because the acceleration is real], we obtain 



a p {u>)* 



(z)t=o 



zo — ft— Re .7-" 
op 



d'Kpc i 



p=p,p±=o 

^a p ((j)* J^OpM • (92) 



2tt 



The forward-scattering amplitude T arises in the one- 
loop diagram without photon emission. Therefore, the 
second term in Eq. (|92|l comes from the quantum cor- 
rection to the potential rather than from the reaction to 
radiation. (We will discuss this term in the next section.) 
Thus, the position shift due to radiation can be given as 



r I Quantum 



ie 2 f +co du 
Girpc 3 2?r 

dttP r {t) , 



d 

a p (u)* — a p (ui) 



P J- 



(93) 



where P r {t) given by Eq. (|16fl is the power radiated. 

One can readily see that Eq. (|93l) does not agree with 
Eq. (|20|) derived in the Lorentz-Dirac theory. However, 
these equations would be identical if we removed the first 
term proportional to \og(v f/vi) from Eq. I|2U|) . This 
term can be traced back to the term proportional to 
vv in Eq. I|1U|) which has arisen due to the substitution 
vv — d(vv)/dt — v 2 . Therefore, Eq. (|93|1 is reproduced in 
the classical theory by modifying the energy-conservation 
equation as 



d_ 



V(z] 



6-7TC 3 



-Pr(t) 



(94) 



This formula can be interpreted as stating that the en- 
ergy of the charged particle is lost through radiation at 
each time. In other words, the ft — > limit of scalar QED 
at order e 2 can be reproduced by classical electrodynam- 
ics if momentum conservation is disregarded and energy 
conservation is used at each moment of time, as far as 
the position shift is concerned. This seems reasonable 
because in the quantum theory the energy is conserved 
in the rest frame determined by the potential V{z) but 
the momentum is not conserved since the potential is 
z-dependent. 



The imaginary part IvaT of the forward-scattering 
amplitude is related to the photon emission probabil- 
ity by unitarity, as we mentioned before. The real part 
Re T represents the quantum correction to the potential 
(which may be nonlocal). To illustrate this fact let us 
compute the forward-scattering amplitude due to a slight 
change in the potential in the form V(z) — > V(z)+5V(z), 
with SV(z) regarded as perturbation. We find 



T 



2mc 



d 3 P 



/+oo 
dz5V(z) 
-oo 



(95) 



to lowest order in V(z) and 5V{z). The position shift 
due to this change in the potential is 



op 



V 



dzSV(z). 



(96) 



It can readily be verified that the position shift is given 
by the same formula in the classical theory. We have 
not been able to find Re J 7 to order h° in our model but 
will show that the quantum correction to the potential 
is of order h . This will imply that the electromagnetic 
correction to the potential at order e 2 cannot be obtained 
in the classical theory. 

Since the potential is part of the mass term in our 
model, its quantum correction is related to the mass 
renormalization. If we let V(z) = in the mass func- 
tion given by 

[m 2 c 2 + 2mV(z)]/h 2 = [m(z)] 2 c 2 /h 2 , (97) 

its quantum correction can be computed in the dimen- 
sional regularization with D = 4 — 2e as 



5m 2 



2 ( 3 3 + 7 31 

167T 2 fi,C \6 47T/U 2 



(98) 



where \i is the renormalization scale and 7 is Euler's con- 
stant. The mass function is corrected according to this 
formula for large and positive z for which V(z) = 0. 
For large and negative z, the squared mass is m 2 + 
2mV- 00 /c 2 . Therefore the quantum correction to the 
mass term is obtained by replacing m by this value in 
Eq. I|98|) . This implies that the difference between the 
quantum corrections to the potential at z = ±00 is of 
order H . Hence, this correction cannot be obtained in 
classical electrodynamics [2l| . 

If the potential is slowly varying so that the WKB 
approximation is valid, it is reasonable to expect that the 
quantum correction to the mass term at leading order in 
ft is obtained by replacing m 2 by [m(z)] 2 in Eq. Ij98(l for 
all z. Then, 



8[m{z)f 



i(z)f 



I6n 2 hc 



e 



7 - 3 In ■ 



i(z)} 2 



Ann 2 



(99) 
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The quantum correction to the potential can be expressed 
by using this formula and Eq. (|98|l as 



V(z) 



3e 2 



m 
[m{z)\ 



32ir 2 hc 



In- 



i{z)f 



(100) 



to lowest order in e 2 . Since the acceleration is —V'{z)/m, 
this formula shows that the correction to the acceleration 
is ultraviolet finite at leading order in h. 



where a is the typical acceleration and where t a is the 
duration of the acceleration p3 |. Then, with p = mv, 



Az 



— >- 

Ap ~ a \at a 



%a ^ > %a ■ 



(103) 



where z a = vt a is the length of the interval in z where 
the acceleration occurs. Thus, the wave packet would be 
much wider than the region of acceleration if the position 
shift were to be larger, in theory, than the width of the 
wave packet at a later time. 



VI. CONCLUSION 

In this paper we studied a wave packet of a charged 
scalar particle moving in the z-direction and accelerated 
by a potential which depends only on z. Our main con- 
clusion is that, in the limit fi->0, the position shift due 
to the reaction to radiation agrees not with the formula 
in the Lorentz-Dirac theory but with that obtained by 
assuming that the kinetic energy is lost to radiation ac- 
cording to the Larmor formula at each moment of time, 
at leading order in e 2 and in the nonrelativistic approx- 
imation. It will be interesting to study a more general 
system and see if the result obtained here can be gener- 
alized. 

Finally, let us discuss the magnitude of the position 
shift in comparison with the width of the wave packet. 
As we have seen in the introduction, the position shift 
may be much smaller than indicated by the fact that it is 
of order h°. First we note that the first term in Eq. Q17[l. 
which represents the discrepancy between the Lorentz- 
Dirac theory and the quantum theory, is in fact much 
smaller than the Compton wavelength of the particle: 



e 2 ^ Vf _ 2ah 
6ttitic 3 ^ ^ Vj 3mc 



miog^ 

V C / V 



(101) 



Thus, it is much smaller than the width of any nonrel- 
ativistic wave packet. Hence, this term would be prac- 
tically unobservable even if it was present. The second 
term, which is the position shift in the quantum theory, 
grows linearly as a function of the time |Tq| after the 
acceleration has taken place. For large |To| this term is 
approximately — E r \To\/p, where E r is the classical en- 
ergy of the radiation. If the width Az of the wave packet 
is narrowest with Az ~ h/ Ap when the particle is going 
through acceleration p^. then after time |To| its square 
behaves like (Az) 2 ~ (h/Ap) 2 + (Ap/m) 2 \T \ 2 . The sec- 
ond term dominates for large enough |Tb|. Therefore, in 
theory, one can make the position shift larger than the 
width of the wave packet by making \Tq\ large and letting 
Ap satisfy 



Ap < 



mE r ahma 2 t„ 



P 



pc* 



(102) 



VII. ERRATUM 

In finding the h — > limit of the scattering ampli- 
tude [Eq. we erroneously discarded surface terms. 
We correct this error here. This correction changes our 
conclusion completely: the position shift in the Lorentz- 
Dirac theory agrees with that in quantum field theory in 
the h — > limit. 

Since the function K(z) in Eq. I|43|) can be approxi- 
mated in the limit h — > by ui/v p (z) [see Eq. (|55() ] we 
have 

/+OO / rZ \ 

dzexp I i / u>/vp(z')dz' 
-co V Jo / 

r + co 

= 2p / v p (t)e iut . (104) 



(Note that P — > p in the h — > limit.) Since this integral 
is ill-defined, we insert a damping factor x(t) which takes 
the value 1 while a p (t) = v' p (t) ^ and goes to zero 
smoothly as t — > ±oo. Thus, 

/+oo 
v p (t) X (t)e^ . (105) 
-oc 

The quantity a p {uj) in Eq. (|89|) needs to be replaced by 



+ OC 



v p (t) X (t)e^ = fipM 



v p (t) x '{ty 



Then the contribution of radiation reaction to the posi- 
tion shift is 



\ rad 



e 2 f + °° dw 
67TC 3 ./ „ 2tt 



X 



a p (—u) + 



v p (t') X '(t')e- iut 'dt' 



+ °° dv (t) 



dp 



6irc 3 



+ 30 



dt 



1 M dv p (t) d 2 

2 v ^^p-dM t)] 



(106) 
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The quantity dv p (t)/dp can be found as follows. Note 
that 



t = 



^p 1 - 2mV{C) 



and hence 
dt = 



-dZn 



V P (t) P rn yj [ Vp ( T )f 
This can be re-arranged as 

* dr 



dp . 



V f 

dz p = v p (t)dt H — -v p (t) 



/o [v p (t)\ 

By using equality of mixed partial derivatives, 



we obtain 



( dz p\ 




( dz p \ 


\dt) 


dt 


\d P ) 



d 
dp 



m 



1 f* dr 

MT) +ap{t) J [v p (T)f 



By substituting this in Eq. (|106fl we find 

f +00 



\ rad 
)t=0 



e 2 Vf 
1 v p (t) 



dt 



z p (t)y 



dr 



[v p (t)Y 



ld_ 

2dt 



a p (t) 



dr 



[v p (t)Y 



e 2 Vf fj vj_ 
Qnmc 3 1 Vi 



dr 



[a p (t)fdt 



(107) 



(108) 



dp. (109) 



(110) 



(111) 



(112) 



where we have used the assumption that a p (t) = if t > 
0. This is exactly the formula obtained by integrating Eq. 
(|12fl found using the Lorentz-Dirac theory from t = — 00 
to t = 0. Notice that we did not use the assumption that 
2m\V(z)\ <p 2 . 
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APPENDIX A: CALCULATIONS WITHOUT 
CHOOSING PHYSICAL POLARIZATION 
STATES 

If we do not use physical polarization states for photons 
and stick to the Fock space with indefinite metric for the 
Feynman gauge, the one-photon final state is not given 
by Eq. (25 but by 



l/,l7) 



d 3 k 



d 3 P 



d 3 p 



2huj{2Tif J 2P (2Trh) 3 J V"2po(27rfi) 3 
x/(p)^(k;P;p)&t( k )At(P)|0), (Al) 



where 

^(k;P;p) = -iehc 2 J dt d 3 x$ P (i,x)* d '„ * P (t,x). 

(A2) 

One may assume that Pj^ and are negligibly small as 
before. Let us define G^{k z ,P,p) by 

^(k,P,p) 

= ec 2 Gp{k Zl P lP ) 

x(2Trh) 3 6 2 {P ± + hk ± - p±)S(P + hu-po) . 

(A3) 

It is clear that G x and G y are negligibly small and that 
Go and G z are independent of the direction of k^. We 
have 



G z (k z ,P,p) = G(k z ,P,p), 



(A4) 



where G(k z ,P,p) is defined by Eq. g2J. The 0- 
component Go{k z ,P,p) is 



G (k z ,P,p) 



Po + Po 



dz 



[k p (z)k p (z)] 1/2 



x exp 



% \ K{z')dz' 



(A5) 



where K(z) is defined by Eq. Q44|l . Using a formula sim- 
ilar to Eq. I]45p. we arrive at 



G (fc z ,P,p) 



-i\J Pp ^° — — - I il: 



dz 



{[n P (z) Kp (z)} 1/2 K(z)} 1 



x exp 



i / K{z')dz' 



(A6) 



At lowest order in h we can let P = p, Pq — po, np{z) = 
K p (z) = mvp(z) and use the approximation l|54(l . Thus, 
we find 



G (k z ,P,p) 



2ipp 



dz- 



k z v'{z) 



c J -00 m [lj — v p (z)k z ] 
exp |i J [lo/v p (z') — k z ] dz'X 



(A7) 
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Then, by letting po — mc 2 and k z = on the right-hand where 
side [because w 3> v p (z)k z ] and comparing the result with 
Eq. lf5fi|l we find 



G (k z ,P, P ) = - — G(k z ,P,p). (A8) 
to 

Thus, the only difference between the calcula- 
tions here and those in Sees. II I II and IIVI is that 
the factor \€^' z (k)G(k z , P,p)\ 2 is replaced here by 
-G„(k„P,p)Gi*{k t ,P,p) = (1 - fc 2 c 2 / w 2 )|G(^,P,p)|' 2 
in the k-integral for the final results. For a given ui the 
average of 1 — k 2 c 2 /uj 2 over the solid angle is 2/3. Since 
the average of |e^^ z (k)| 2 is also 2/3, we arrive at the same 
transition probability l|77(l and position shift Ij93(l . 



APPENDIX B: A DERIVATION OF THE 
POSITION EXPECTATION VALUE 

In this Appendix we derive Eq. I|80|l which was used to 
find the expected value of the z-coordinate of the particle. 
A systematic method to obtain a wave function for a one- 
particle state in the scalar field theory has been given 
by Feshbach and Villars 0, |2(J ■ We will not use this 
method explicitly, but the derivation we describe here is 
based on the idea spelled out in Ref. [lfij . 

Recall that the charge density operator is 



(Bl) 



&(t,x) = — <0|V(t,x)|i) 
d 3 V 



: /(p)* p (t,x). (B4) 



Then it follows from orthonormality (|32|l of $ p (i,x) and 
the normalization condition Ij67() for /(p) that 



J d 3 xpi(t,x) = 1. 



(B5) 



Thus, the function pi(t, x) is the normalized probabil- 
ity density function for the position of the particle. By 
substituting Eq. ljB4|) in Eq. I|B3|) we obtain 



Pi(*,x) 




1 f d 3 p 

2 J {2irh) 3 J {2nh) 3 \\jp Q 
xf(p'yf(p)%>(t,xy%(t,x), (B6) 

where po is defined by Eq. i|26|) with p' Q similarly defined. 
The expected value of the z-coordinate at t = is 



(B7) 



Let us define 



By using the assumption that the wave packet is located 
far into the region where V(z) = at i = we may 
approximate the mode function <J> p (0,x) at t = in Eq. 
(151) as 



pi(t,x) = (i\p(t,x)\i) 



(B2) 



where the state |z) is defined by Eq. (|66|l . Since this 
state is a one-particle state, the charge density pi(t,x) 
coincides with the probability density function for the po- 
sition of the particle at each time up to a normalization 
factor (which will turn out to be one). By a straightfor- 
ward calculation we find 



Pi(t,x) = ihcj>i{t,x.)* d t <f>i(t,x) , 



(B3) 



$ p (0,x) w exp 



(pz + p ± ■ Xj 



With this approximation we can write 



d 



z$ p (0,x) « -ih— $ p (0,x). 



(B8) 



(B9) 



By substituting this formula in Eq. I|B7(1 and integrating 
by parts we find 



ih 



= ih 



d 3 p f d 3 p' d 



2 ./ (2nh) 3 J (2Trh) 3 dp 
d 3 p f d 3 p' 



^ + J-)/(p')*/(p) 
Po V Po 



y"d 3 x$ p ,(0,x)*$ p (0,x) 



(2iTh) 3 J (2irh) 3 
_ 1 Po 



2 ^- 7 _ w 5)/ip'.7(p) + /(p') -Hp) 



d_ 

dp J 



(2ttH) 3 S 3 (p - p') 



= ih 



d 3 r> d 



(2nh) 3J ^' dp J 

We immediately obtain Eq. I|8U|) from this formula by integration by parts. 
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